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Three-dimensional (3D) photonic crystals with a 3D photonic bandgap play a fundamental role 
in cavity quantum electrodynamics (QED), especially in phenomena where the local density of 
optical states is essential. We first review the current status of the fabrication of 3D photonic 
crystals with a bandgap at optical frequencies, corresponding to wavelengths below 2500 nm. We 
discuss the main implications of 3D bandgaps for cavity QED, in particular spontaneous emission 
inhibition of emitters embedded in a 3D bandgap crystal. Moreover, we review progress in enhanced 
emission of emitters placed in a photonic bandgap cavity, thresholdless laser action in a miniature 
photonic crystal cavity, breaking of the weak-couping limit of cavity QED. Finally, we discuss 
several exciting applications of 3D photonic band gap crystals, namely the shielding of decoherence 
for quantum information science, the manipulation of multiple coupled emitters including resonant 
energy transfer, lighting, and possible spin-off to 3D nanofabrication for future high-end computing. 


INTRODUCTION 


The propagation of light in periodically ordered pho¬ 
tonic crystals bears a strong analogy to the wave prop¬ 
agation of conduction electrons in atomic crystals. The 
development of Bloch modes and the dispersion are de¬ 
termined by the interference of waves that are diffracted 
by lattice planes [H, i]. The periodicity commensurate 
with the wavelength (a ~ A/2) gives rise to Bragg diffrac¬ 
tion that is associated with frequency windows for which 
waves are forbidden to propagate in a particular direc¬ 
tion. Such stop gaps have long been known to arise for 
electromagnetic waves, notably for X-rays in atomic crys¬ 
tals 0. At optical frequencies, the stop gaps that occur 
in one dimensional periodic structures, known as Bragg 
stacks, are widely used in optical laboratories as broad¬ 
band highly reflecting dielectric mirrors Q. 

The distinguishing feature of three-dimensional (3D) 
photonic crystals is that a common stop gap can be 
achieved for all directions and for all polarizations si¬ 
multaneously: the widely pursued 3D photonic bandgap. 
Historically a 3D bandgap has never been considered for 
the propagation of X-rays in periodic media, since a gap 
requires a high refractive index contrast of order unity, 
whereas the refractive index varies by less then 10“"^ in 
the X-ray range. At frequencies inside the 3D photonic 
bandgap, the density of optical states (DOS) completely 
vanishes. As the density of states can also be interpreted 
as the density of vacuum fluctuations a 3D bandgap 
thus serves as an effective shield to these fluctuations. 
The total absence of optical modes in a photonic bandgap 
has implications beyond classical optics that break the 
analogy between the behavior of light in photonic crys¬ 
tals and the behavior of electrons in atomic crystals. 

3D photonic bandgap crystals play an important role 


in cavity quantum electrodynamics (cQED) @,0 where 
they offer at least five prospects for new physics Q. 
Firstly, probably the most eagerly pursued phenomenon 
is the complete inhibition of spontaneous emission: an 
excited quantum emitter - such as an atom, molecule, or 
quantum dot - embedded in a crystal with its transition 
frequency tuned to within the 3D bandgap remains for¬ 
ever excited since it cannot decay to the ground state 
by emitting a photon. Any interaction mediated by vac¬ 
uum fluctuations is affected by their suppression in the 
3D bandgap 0. Therefore, a crystal with a 3D pho¬ 
tonic bandgap not only inhibits spontaneous emission - 
including a shift of the emitter’s frequency known as the 
Lamb shift [13 - it will also modify the spectrum of black- 
body radiation ll|, it will affect resonant dipole-dipole 
interactions including the van der Waals and Casimir 
forces [13 113 , and the well-known Forster resonant en¬ 
ergy transfer that is prominent in biology and chem¬ 
istry dill. 


A finite inhibition of spontaneous emission is also fea¬ 
sible in other optical materials such as microcavities or 
nanowires, if the emitters are carefully positioned in a 
tiny volume [mil- In contrast, in photonic bandgap 
crystals the extent over which emitters are controlled is 
only limited by the crystal’s volume. Moreover, the com¬ 
plete and radical suppression of electromagnetic modes 
in the bandgap is unique to photonic crystals, and is not 
found in other optical materials that appear to exclude 
all light in a particular frequency range. For instance, if 
we imagine a metal-coated box, then an emitter inside 
such a box is not seen from the outside, and the emission 
might be perceived to be inhibited. Nevertheless there 
are many optical states in the box in which photons can 
be emitted, before they are ultimately absorbed in the 
metal walls. In a 3D photonic bandgap, however, there 
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are simply no electromagnetic states available hence an 
excited emitter cannot emit a photon at all and remains 
forever in the excited state. 

Secondly, once a bandgap is achieved, the cQED 
physics becomes even richer by introducing point defects. 
Such defects locally break the crystal’s symmetry, which 
results in the appearance of isolated electromagnetic res¬ 
onances in the bandgap. At these resonances the field is 
spatially localized within a tiny nanoscale volume Vcav 
that is typically less than a wavelength cubed and thus 
less than a cubic micron. In other words, a point de¬ 
fect acts as a tiny cavity that is shielded in all three 
dimensions from the vacuum by the surrounding crys¬ 
tal 17, Hence such a photonic bandgap cavity is 
called a ” nanobox for light”. Since the density of states 
in a nanobox is proportional to IjVcav and thus greatly 
enhanced by the tiny volume, an embedded emitter ex¬ 
periences a greatly enhanced emission rate, also known 
as Purcell effect [l9j |. 

A third reason why 3D photonic bandgaps are rele¬ 
vant to solid state cQED occurs when a gain medium is 
introduced in a nanobox. Such a nanobox with gain of¬ 
fers the promise of a thresholdless laser Q: Since only 
one resonance exists in a photonic bandgap cavity, and 
the vacuum is shielded, there is no competing sponta¬ 
neous emission into modes other than the lasing mode, 
and the laser thus immediately switches on [i^. More¬ 
over, since 3D photonic bandgap crystals are typically 
semiconductor devices, the on-chip integration of such a 
thresholdless laser is readily foreseeable. 

Eourth, an important research theme in cQED is the 
breaking of the weak-coupling approximation. There are 
in essence two ways to break this limit: The most well- 
known approach consists of embedding a two-level quan¬ 
tum emitter in a high-finesse cavity, and tuning the emit¬ 
ter frequency ujeg to the cavity resonance frequency io^av ■ 
When a quantum of energy is exchanged between the 
emitter and the cavity at a rate D/j - called the vacuum 
Rabi frequency - that exceeds leakage rates such as the 
spontaneous emission rate or the cavity escape rate, the 
emitter and the cavity resonance are in a coherent su¬ 
perposition, called QED strong coupling While the 
achievement of this limit has been realized with pillar 
cavities 21|^ing cavities (2^, and 2D photonic crys¬ 
tal cavities [23|, the realization in nanoboxes is still out¬ 
standing. There are several alternative ways to break the 
weak-coupling limit. One approach is to operate close to 
a van Hove singularity where the density of states has a 
cusp 0. A second approach to break the weak-coupling 
limit consists of rapidly modulating the ” bath” that sur¬ 
rounds a two-level quantum emitter [2^, using ultrafast 
all-optical switching methods [i^ . 

Eifth, in quantum physics there is an active interest 
in decoherence, that is, the loss of coherence between 
the components of a system that is in a quantum su¬ 
perposition [ 2 ^. A consequence of decoherence is that a 


quantum system irreversibly reverts to revealing classical 
behavior, which is undesirable for applications in notably 
quantum information processing [27|. In case the quan¬ 
tum systems are optical emitters, cQED is the relevant 
field. An important component of decoherence in cQED 
corresponds to the escape or emission of photons, or by 
absorption of the photons by the environment 28| . Hence 
the shielding of vacuum fluctuations by a 3D photonic 
bandgap offers opportunities to make optical quantum 
systems robust to decoherence. 

In this review we summarize recent progress on the five 
subjects in cQED listed above, with emphasis on experi¬ 
mental work. In addition, we provide an overview of the 
current status regarding the fabrication of 3D photonic 
bandgap crystals. For the scope of this review, we limit 
ourselves to i) classes of photonic crystals with demon¬ 
strated 3D photonic bandgaps, and ii) light frequencies in 
the optical regime, chosen to correspond to wavelengths 
A < 2500 nm (or frequencies u:/{2tt) > 10^“^ s“^) as these 
are accessible by conventional optics available in many 
laboratories all over the world. Moreover, at these fre¬ 
quencies stimulated and spontaneous emission rates are 
highly significant. 


THEORY 


Spontaneous emission control 


In the weak-coupling approximation of quantum elec¬ 
trodynamics (QED), the radiative rate of spontaneous 
emission ^rad of a two-level dmole quantum emitter is 
given by Fermi’s’’golden rule” [29j. It is well-known that 
spontaneous emission is not an immutable property of an 
emitter - such as an excited atom, molecule or quantum 
dot - but it is also influenced by the emitter’s nearby en¬ 
vironment [IS 113. The influence of any non-dissipative 
environment is described by the local density of opti¬ 
cal states (LDOS) N{oj,r,ed) that counts the number of 
photon modes available for emission weighted with their 
amplitude squared, and that is interpreted as the density 
of vacuum fluctuations [3l|. Control of QED properties 
by means of confined light - that is expressed by a mod- 
ihed LDOS - is generally considered to be the realm of 
cavity QED [3 • The radiative rate for the transition from 
the excited state IQ to ground state \g) is conveniently 
expressed as IJ, l32l| 


^radi^eg — 


^eg 


heo 


N{uj = ujeg,r,ed), (1) 


with LOeg the emission frequency, r the position of the 
emitter, the dipole orientation, d the modulus of the 
transition dipole moment matrix element, and h Planck’s 
constant. It is instructive to briefly summarize the ap¬ 
proximations and assumptions used to derive the time- 
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evolution of the emitter’s exited state and thus Fermi’s 
” golden rule” Eq. ([T]) : 


a. The electric-dipole approximation is applied to the 
emitter to evaluate the electric-field operator. 

b. The perturbation expansion is used to describe the 
time evolution of the emitter-field system. 

c. The rotating-wave approximation is used to neglect 
rapidly-changing counter-rotating field terms. 

d. The dielectric function e(r) is taken to be real to 
properly define modes and hence the LDOS (see 
Eq. [3] below.) 


e. The emitter-field interaction is taken to be a 
Markovian process, which assumes that if a photon 
is emitted the memory of this event is lost practi¬ 
cally instantaneously by the quantum systerm thus 
corresponding to a vanishing coherence time |5l.l24|. 


From equation (HD it is apparent that one can control 
the emission rate by means of the LDOS. The prefac¬ 
tor contains intrinsic emitter properties namely the tran¬ 
sition dipole moment d and transition frequency ujeg- 
Equation (HD reveals the well-known fact that the emis¬ 
sion rate depends on the frequency and the position of the 
emitter. While many efforts are directed at controlling 
spontaneous emission by various classes of nanophotonic 
systems - cavities, antennae, plasmonic systems, random 
media - arguably the most radical change in emission oc¬ 
curs in a 3D photonic bandgap where the LDOS vanishes, 
see FigurelH Therefore, an emitter in a 3D bandgap that 
is in the excited state is radically forbidden to decay by 
emitting a photon. In other words, while the excited 
state |e) is usually unstable against decay to the ground 
state \g) under the emission of a photon, the excited state 
is stabilized in a 3D photonic bandgap! 


The local density of optical states 

The local density of optical states (LDOS) N{u},r,ed) 
is defined as: 

N{uj,r,ed) = —(ej • Im(G(a;; r, r)) • e^), (2) 

'Trn^ 



Reduced frequency {a/x) 


FIG. 1. Local density of optical states (LDOS) versus fre¬ 
quency in a 3D photonic bandgap crystal. The LDOS was 
calculated for an infinite inverse opal with an e = 11.9 back¬ 
bone (red connected circles), relevant for silicon. The LDOS 
was averaged over the unit cell such that it represents the total 
density of states. The blue hatched bar indicates the range 
where the LDOS vanishes: the 3D photonic bandgap. The 
black dashed curve is the LDOS of a homogeneous medium 
with an effective refractive index similar to that of the crystal. 
The frequency is reduced with the lattice constant a. Data 
from Ref. 


is more tractable to calculate the LDOS by a summation 
over all Bloch modes as 3l|: 


1 

N{uj,r,ed) = . .3 , . ^ / dk(5(w-Wk,p)|ed-Ak,p(r)p, 

fzTTj efrj ^ J -00 

( 3 ) 

with Ak,p(r) a field mode with wavevector k and polar¬ 
ization state p = 1,2. For cavity QED it is important 
to note that the field modes are projected on the ori¬ 
entation of the transition dipole moment e^. Since this 
expression for the LDOS employs the plane-wave expan¬ 
sion [i^, the LDOS in Eq. ([3]) pertains to an infinitely 
extended crystal L —>■ 00 . In the limit of a homogeneous 
dielectric medium with a spatially independent refractive 
index n = -y/e, the LDOS is equal to 


N{u}) 


(37r)2c3’ 


( 4 ) 


with G(w; r, r) the Green dyadic 0 that depends on fre¬ 
quency and position. The LDOS is a classical quantity, 
which can be appreciated from the absence of Planck’s 
constant from Eq. ([2]). In case of dissipative optical me¬ 
dia with complex e, the Green dyadic is well defined, and 
the imaginary part of the Green dyadic describes the to¬ 
tal decay rate, Le., the sum of the radiative decay rate 
and the non-radiative quenching rate induced by the dis¬ 
sipative environment. In dielectric photonic crystals, it 


which illustrates the well-known dependence on the fre¬ 
quency squared. 

Figured] shows the LDOS calculated for a 3D photonic 
bandgap crystal. At low frequencies the LDOS increases 
quadratically, which illustrates the feature that if the 
wavelength is much greater than the lattice parameter 
(A » a), the crystal effectively behaves as a homogeneous 
medium, c/. Eq|3| With increasing frequency modula¬ 
tions appear in the LDOS, as well as characteristic cusps 
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called van Hove singularities 0. The blue hatched bar 
highlights the range where the LDOS completely van¬ 
ishes: the 3D photonic bandgap. The complete vanish¬ 
ing of the LDOS over a finite bandwidth, irrespective of 
position r or dipole orientation e^, is a unique feature of 
3D photonic crystal. Whereas other systems such as 2D 
slab photonic crystals Q or nanowires 161 have deep 


pseudogaps and concomitant strong emission inhibition, 
the LDOS never completely vanishes. 

At the edge of the photonic band gap the LDOS van¬ 
ishes as ^J\uj — LOedge \ as notably pointed out in Ref. . 
Hence the edge has a cusp and is thus also a van Hove 
singularity. For many years, the LDOS was considered 
to diverge at the edge of the gap, which led to many in¬ 
tricate QED predictions [36|. It appears, however, that 


such a diverging LDOS is typical for ID systems [10| , but 
not for 3D crystals. 


Quantum efficiency of the emitters and degree of 
cQED control 


In order to choose a suitable QED experiment with 
emitters in a photonic bandgap crystal or to correctly 
interpret results, it appears to be essential to consider 
at least one property of the emitters, namely the emis¬ 
sion quantum efficiency 37|, [3^ . The emission quantum 


efficiency rj of one emitter is defined as the ratio of the 
radiative rate to the sum of the radiative and the non- 
radiative rate ^nrad- 


Vi^eg) = 


^radi^eg) _ '~)rad i^eg) 

'Jradi^eg'} ^nrad ^toti^eg') 


(5) 


For the experimentally relevant case of an inhomoge- 
neously broadened ensemble of emitters, the quantum ef¬ 
ficiency becomes distributed; this situation is extensively 
discussed in Ref. [s^. 

High-efficiency quantum emitters. If the quantum ef¬ 
ficiency of the emitters is high (77 —>■ 1 ), the experi¬ 
mental method of choice is time-resolved emission. In 
this method, the decay of the population density of ex¬ 
cited emitters N^xcit) is probed by recording a decay 
curve g{t). Here g(t) is the total decay, i.e., the sum 
of the radiative decay events f(t) and non-radiative de¬ 
cay events {g — f){t). We are primarily interested in 
time-resolved emission measurements, which are gener¬ 
ally recorded by the well-known time-correlated-single- 
photon-counting method . The resulting decay curve 
f{t) is the distribution of arrival times of single photons 
after exciting the emitters with a short laser pulse at 
time t = 0 , averaged over many excitation-detection cy¬ 
cles. Such a histogram is the probability density of emis¬ 
sion which is modeled with a probability density func¬ 
tion The emission by the excited emitters at time 


t is described by a reliability function or cumulative dis¬ 
tribution function equal to ) — !■ relia¬ 


bility function tends to I in the limit t —>■ 00 and to zero 
in the limit t ^ 0. The relation between the fraction of 
excited emitters and the decay curve, in other words, be¬ 
tween the reliability function and the probability density 
function is 


f g{t')dt' = 1 - 
^0 


Nexcjt) 

NexM' 


( 6 ) 


Physically equation m means that the decrease of the 
density of excited emitters at time t is equal to the in¬ 
tegral of all prior decay events [s^. Equivalently, the 
total decay g{t) is proportional to the time derivative 
of the fraction of excited emitters. In many reports the 
distinction between the reliability function and the prob¬ 
ability density function is neglected; the intensity of the 
decay curve g{t) is assumed to be directly proportional 
to the density of excited emitters Nf.xc(,t). This propor¬ 
tionality only holds for single-exponential decay but not 
for the general case of non-single-exponential decay; this 
distinction has important consequences for the interpre¬ 
tation of non-single-exponential decay that often occurs 
in photonic crystals (39|. 

From the slope of the photon arrival time distribu¬ 
tion, we obtain the total decay rate that is equal to 
Itot = Irad + Inrad- Since an efficient emitter has 
Irad ^ Inrad, 8t time-resolved experiment yields to a 
good approximation the radiative rate jrad that one is 
interested in. Since the radiative rate depends on the 
density of states, see Eq. o, it is generally independent 
of the direction of emission. While there are special situ¬ 
ations that the rate may depend on detection angle (see 
Ref. [42|), one should beware of possible artifacts in this 
respect. 

The average arrival time of the emitted photons or the 
average excited-state lifetime is given by the first moment 
of the time-resolved emission curve f{t) [ 3 ^. This result 
confirms that the dynamics of the population of the ex¬ 
cited state Nexc{t) is controlled by embedding emitters 
in a photonic band gap crystal. 

Low-efficiency quantum emitters. If the quantum effi¬ 
ciency of the emitters is low (77 <C 1), the experimental 
method of choice is the continuous-wave (cw) observation 
of the emitted intensity 37,|38|. From rate equations one 
can derive that the emitted intensity Iffiieg) is equal to 


I{wieg) = 'ytot{uJeg)Nexcri{^eg) = Plji^eg) = ^ 

'ytotyOJeg) 

where P is the rate of excitation of the emitters [^ . 
If the quantum efficiency is low then jrad ^ Inrad, 
hence jtot —t Inrad so that ■jtot is independent of the 
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LDOS. Consequently the cw intensity I{oJeg) of the emit¬ 
ters becomes proportional to yrad, hence I{uieg) is then 
an excellent probe of the LDOS. If one studies the cw 
intensity I (to eg) to obtain information on the crystal’s 
LDOS, care must be taken to distinguish spectral fea¬ 
tures from several other effects, such as angular effects 
related to photonic bandstructures, as well as the angle- 
dependent excitation and the angle-dependent collection 
efficiency [ailii- To interpret results in terms of the 
LDOS, the observed intensity spectrum I{ujeg) is normal¬ 
ized to a reference that employs the same light sources 
in the same chemical environment. Hence the same tran¬ 
sition dipole moment and the same quantum efficiency 
pertain, so that equations mum are properly inter¬ 
preted. Second, a reference must have a well-known 
LDOS. Third, a reference is required with an environ¬ 
ment where the emitters experience the same refractive 
index as in the photonic bandgap crystal, to avoid com¬ 
plications arising from Lorentz local field factors g. 
Fourth, since cw intensity is very sensitive to the collected 
solid angle, it is important that both the crystal and the 
reference have similar escape functions for internal light 
escaping to the detection system, or at least very well 
known escape characteristics. It is our experience that 
the best suited reference is a similar photonic crystal in 
the long wavelength limit, for several reasons: the LDOS 
is well-known and proportional to , the average index is 
the same hence Lorentz local field factors cancel, and the 
escape function is well-known and matches a Lambertian 
distribution [s^. Finally, this choice ensures a reference 
made by the same fabrication process, hence the same 
chemistry and thus the same non-radiative rate. Other 
considered reference systems such as a random medium, 
or a bulk material do not share these strengths and are 
therefore less reliable reference systems for the interpre¬ 
tation of the cw intensity. 

In case of low-efficiency quantum emitters, the total de¬ 
cay rate measured in a time-resolved experiment hardly 
depends on the LDOS as it is dominated by non-radiative 
decay events ((ff—/)(t)) [s^l- In this situation the dynam¬ 
ics of the excited-state population is hardly controlled 
by the photonic environment, in contrast to the case of 
efficient emitters whose population Nexc(t) is truly con¬ 
trolled by the photonic band gap environment. 


Beyond weak coupling 


[ht] In the Markov approximation we assume that the 
quasi-continuum of the field modes - the hath - equili¬ 
brates on time scales much faster than the oscillation pe¬ 
riod of the dipole r^. Since the LDOS iV(a;,r,ed) has a 
dimension of per frequency, it denotes a timescale that is 
usually considered to be the correlation time of the bath 
Tb 2M ■ From Eq. U we calculate the correlation time 


for the LDOS at an optical frequency of 10^° s ^ in a 



Scaled time (t.n„ or t.y J 

' R 'rad' 


FIG. 2. Excited state population Ne^c of a two-level emitter 
with time t. In the strong-coupling limit, damped vacuum 
Rabi oscillations appear (red full curve). In the weak-coupling 
limit, the population of the excited state |e) decays exponen¬ 
tially with the spontaneous emission rate 'jrad (black dashed 
curve). Time is scaled with the vacuum Rabi frequency Q.r, 
or the spontaneous emission rate yrad in case of weak cou¬ 
pling. In the strong-coupling case, 7rad is taken to be equal 
to Qr. 


quantization volume V = (27rc/a;)^ to be u = 10“^® s, 
which is much shorter than Td — 10“^® s. The interac¬ 
tion with the fast bath destroys the emitter’s memory 
of the past, which leads to the exponential decay of the 
excited state. The situation ^ Td is also known as 
the weak-coupling limit of cavity QED Bi, and per¬ 
tains to the majority of spontaneous emission studies in 
nanophotonics with photonic crystals. 

Highly interesting quantum physics arises when the 
weak-coupling limit is violated. When the LDOS is 
strongly increased, the relaxation of the bath becomes 
slower. Once the correlation time Tb becomes of the 
order of, or even longer than the dipole oscillation pe¬ 
riod Td, the quantum system retains a memory of ear¬ 
lier times, hence approximation e of the weak-coupling 
limit is violated. [45j| This means that the quantum emit¬ 
ter and the bath are becoming strongly coupled, and the 
evolution of the emitter’s excited state |e) does not show 
exponential decay anymore. In this situation - typically 
studied with a high-Q cavity (2ll - [2^ - the emitter reso¬ 
nance and the cavity resonance hybridize, such that the 
emitters’s excited state population oscillates in time, as 
illustrated in Figure [2l 

In addition to raising the LDOS by a cavity resonance, 
there are other ways to violate the weak-coupling limit. 
One approach is to operate close to a van Hove singular¬ 
ity i- Such a singularity manifests as a cusp in the 
density of states, causing the density of states to be¬ 
come non-analytical. The non-analytic behavior means 
that approximations b and e of the weak-coupling ap- 
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proximations are violated. As a result, a single emit¬ 
ter tuned to a van Hove singularity is predicted to ex¬ 
hibit non-exponential dynamics, including intricate time- 
dependent oscillations of the excited-state population 
that are called ’fractional decay’ Idill^. 

Another approach to violate the weak-coupling limit is 
rapidly modulate the bath. Hence approximations b and 
c of the weak-coupling limit are violated. Rapid modula¬ 
tion of the bath can be achieved by the ultrafast switching 
of the optical properties of a 3D photonic bandgap crystal 
on ultrafast times scales, using all-optical methods (25l |. 


ULTIMATE TOOLS FOR 3D PHOTONIC 
BANDGAP CAVITY QED 

Requirements for 3D photonic bandgap crystals 


There are four main requirements to periodically or¬ 
dered nanostructures in order to successfully function as 
3D photonic bandgap crystals. Firstly, the photonic in¬ 
teraction strength S between the light and the nanostruc¬ 
tures should be elevated. In the introductory chapter, we 
have seen that S is defined as the polarizability per vol¬ 
ume. The photonic interaction strength equals to a very 
good approximation of the relevant frequency bandwidth 
Auj/lo of a stop gap: Auj/uj = S. For photonic crystals 
composed of spherical scatterers, S can be expressed an¬ 


alytically as [48|, It 


^'KOL 6 — 1 


( 8 ) 


where e is the relative dielectric function of the spheres, 
equal to e = vr? with m the ratio of the refractive in¬ 
dices of the spheres and of the surrounding medium: 
m = nijn 2 . Thus a high refractive index contrast m 
is highly desired, which dictates the choice of the con¬ 
stituent materials. Air with index n = 1.0 is a convenient 
low index material, and semiconductors such as silicon or 
GaAs are often the high index material of choice, with 
an index n ~ 3.5. [H^ In Eq. ([8]), g{K.R) is the form fac¬ 
tor of the scatterers (here: spheres in the Rayleigh-Gans 
limit) as a function of scatterer radius R and modulus of 
the diffraction vector AT = |K| = |kout ~ kinj. The form 
factor dictates both the shape of the scatterers, and the 
optimal filling fraction of the scatterers in the crystal [i^ . 

The second main requirement is that the optical ab¬ 
sorption of the constituent materials should be as small 
as possible. This is borne out by the fact that a photonic 
bandgap is a multiple scattering phenomenon, hence at 
every chance that photons are absorbed, it’s ’’game over”. 
The role of absorption has been studied via the imag¬ 
inary part of the frequency w" (blj. Near the edge of 
the gap the square root singularity disappears for the 
weakest absorption considered. Within the bandgap the 
density of states becomes non-zero and proportional to 


Lo"jujgap- Hence for cavity QED in a 3D bandgap, ab¬ 
sorption of light should be reduced as much as possible. 
This requirement limits the use of semiconductors to fre¬ 
quencies below the electronic bandgap. Thus photonic 
crystals made from silicon or GaAs are limited to wave¬ 
lengths longer than 1100 nm or 870 nm. Indeed, in the 
early days of semiconductor research, an empirical rela¬ 
tion between refractive index and electronic bandgap Eg 
was established called Moss’ rule: n^Eg = 77 In the 
visible range, Ti 02 is a versatile high-index material with 
n ~ 2.7, depending on its atomic crystal structure. Thus 
the choice for a high index material limits the bandgap 
width (via the maximum index) and the operating fre¬ 
quency (via the electronic gap). 

The third main requirement is the photonic crystal 
topology. In Ref. it was found that for electromag¬ 
netic waves the network topology is more favorable for 
the appearance of gaps than the Cermet topology with 
isolated scatterers. In the network topology, both low 
and the high-refractive index materials form a continu¬ 
ous network. This requirement strongly limits the classes 
of relevant crystal structures, as for instance the impor¬ 
tant class of isolated scatterers in a medium - typical for 
suspension of colloidal nanoparticles - do not meet this 
requirement. 

The fourth main requirement to photonic crystal struc¬ 
tures is that the unavoidable statistical variations of the 
crystal dimensions - in other words, random disorder - 
must be as small as possible. It is convenient to express 
random disorder in terms of relative variations of the 
lattice parameter 5a/a, and of size variations of the con¬ 
stituent building blocks, such as the relative variation of 
sphere radii 6R/R. Random structural variations cause 
light to be scattered leading to the extinction of coherent 
beams 5^. In case of a high photonic strength typical 
of bandgap formation, the extinction length is limited 
to about 200 lattice spacings in case of small structural 
variations of 1%, typical of a high-quality structure. We 
will see below that significant photonic bandgap effects 
have recently been observed on excited-state population, 
from which it is concluded that state-of-the-art structures 
have sufficiently high order for cavity QED. Therefore, 
we note that a 3D photonic bandgap seems to be more 
robust to disorder than slow light phenomena, which are 
also considered in the context of strong light-matter in¬ 


M|. 


teraction 

A large variety of methods has been proposed to obtain 
many different types of 3D photonic crystals. Examples 
of the fabrication of these structures are described in a 
multitude of reviews. Refs. 18|, offer excellent in¬ 
troductions into the field of photonic crystals where the 
theory of photonic crystals, fabrication methods, and em¬ 
bedding of cavities are extensively discussed. Ref. 
discusses self-assembly as a tool to obtain both ordered 
3D photonic crystals and random structures. Ref. 581 
reviews photonic crystal fabrication by direct laser writ- 
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TABLE I. Overview of three main classes of 3D photonic crys¬ 
tals. The maximum calculated width of the bandgap is given 
with the relevant references. 


on comparing the measured stopbands with calculated 
band structures may be impeded by the fact that the 
calculated geometries potentially differ from real crys¬ 
tals due to inevitable unknowns in the fabrication or the 


Type of 

photonic crystal 


Calculated max p^i^^^^^^^^Jharacterization. Moreover, the complication can arise 
width of the bandgap in Ref. fiel^odes only couple to a specific external polar- 


inverse air-sphere opals {fee) 12% for nsi = 3.45 

Inverse woodpiles (diamond-like) 25% for nsi = 3.6 
Woodpiles (diamond-like) 18% for nsi = 3.6 [59] 


mg. 

This section focuses on fabrication methods for 3D 
photonic crystals with demonstrated bandgaps for light 
with wavelengths up to 2500 nm. A few notable classes 
of photonic bandgap crystals are identified, including 
inverse opals, and diamond-like woodpiles and inverse 
woodpiles, see Table HI In addition 3D photonic crystals 
containing optical cavities are discussed. 


Optical signature of a 3D photonic bandgap 

An important aspect of 3D photonic crystal fabrication 
is their characterization. Scanning electron micrographs 
(SEM) are often used to investigate the outside of the 
resulting structure and obtain a first impression whether 
the fabrication method was successful, see the SEM im¬ 
ages in almost any photonic crystal paper in a materi¬ 
als science journal. To investigate the interior quality of 
the crystals it is possible to “open” the structures after 
their fabrication. For instance, focused ion beam milling 
was employed to (destructively) observe the inner struc¬ 
ture in Refs. [60|,l6]|. A non-destructive characterization 
method for 3D nanostructures is for instance small-angle 
X-ray scattering [^. An experimental demonstration 
that 3D crystals exhibit a photonic bandgap, however, is 
more challenging and requires advanced optical methods 
and a careful analysis. Structures with a 3D photonic 
bandgap have, by definition, the following properties: i) 
stopbands with overlapping frequencies for all wave vec¬ 
tors and polarizations simultaneously, and ii) a, vanishing 
density of states. These properties allow the presence of 
a bandgap to be probed experimentally. 

Experimentally the widths of peaks in reflectivity or of 
troughs in transmission provide a good estimate of the 
width of a stopband. To assess the photonic strength of 
the crystal, the width of the stopbands in experiments is 
often compared to the width of stop gaps from calculated 
bandstructures, see Refs. [63], [6^. Care must be exerted, 
however, as stopbands in reflectivity or transmission may 
also occur due to so-called silent modes. These modes oc¬ 
cur when incident plane waves cannot couple to a field 
mode inside the crystal with a peculiar spatial symme¬ 
try 18l.l65ll66l|. Furthermore, deriving a conclusion based 


“fzation [67| . Thus polarization-resolved experiments are 
necessary to demonstrate that stopbands occur for all po¬ 
larizations [68| . It is also important to exclude spurious 
lioundary effects by verifying that measured stopbands 
reproduce at different locations on the fabricated pho¬ 
tonic crystal. Consequently, for reflectivity or transmis¬ 
sion experiments to serve as a reliable indicator for the 
presence of a bandgap they must be performed i) such 
that all incident angles are probed (27r sr solid angle), 
ii) with resolved polarizations. Hi) on different exposed 
surfaces of the structure, and iv) shown to be position 
independent (69| . 

The main property of photonic bandgap crystals, 
namely a vanishing density of states, allows for a robust 
investigation of the bandgap. However, it is challenging 
to experimentally probe the density of states. The main 
approach is to embed suitable emitters in the crystals 
and measure their emission rates, as originally done for 
inverse opals with a pseudo gap but no 3D bandgap in 
the LDOS In a photonic bandgap the sponta¬ 

neous emission of the emitters is inhibited, thus result¬ 
ing in longer lifetimes of the excited state as described 
by Fermi’s golden rule Eq. ©■ Therefore the lifetime 
is an observable that directly relates to the density of 
states and one can argue that measuring the emission 
rate of embedded emitters is the most suitable experi¬ 
mental method to demonstrate a 3D photonic bandgap. 
Relevant experiments are discussed in the next section 
on cavity quantum electrodynamics. 

When describing how photonic crystals were obtained, 
in many papers the fabrication data are complemented 
with reflectivity and transmission spectra from one or 
a few directions only. In light of the discussion above, 
in a strict sense such results are no experimental proof 
of a bandgap. To give an overview of a wide spectrum 
of available fabrication procedures, we take a pragmatic 
approach and include results in which a bandgap was in¬ 
ferred from such measurements in this Review. Table m 
provides an overview of the fabrication methods and dif¬ 
ferent types of photonic bandgap crystals discussed, in¬ 
cluding relevant references. 


Inverse opals 

In inverse opals, spherical voids are stacked in a 
face-centered-cubic (fee) structure. These air spheres 
are embedded in a backbone material with a high in¬ 
dex of refraction. Typical fabrication procedures for 
these photonic crystals employ template-assisted assem- 
















bly based on templates of close-packed, fcc-ordered col¬ 
loidal spheres, infiltrating the templates with a high re¬ 
fractive index material, and subsequently remo ving the 


template material to obtain air-sphere crystals [71 


For completely infiltrated inverse opals with a high-index 
volume fraction of </) = 26% it was established that the 
refractive index contrast must exceed m = 2.8 in order 
to open a bandgap [t^ . Subsequent calculations revealed 
that intricate incomplete filling of high-index material in 
the form of shells surrounding the spherical voids and 
connecting windows between the spheres yields a max¬ 
imum possible bandgap width of Aui/to = 12% for fee 
silicon inverse opals 4^ 7^. 

Silicon is an excellent backbone material due to its high 
dielectric constant of 11.9 (refractive index n = 3.45) at 
A = 1550nTO and the availability of routine deposition 
methods. Silicon inverse opals were first demonstrated 
in Refs. 7^, lZ6|. In both cases, the template was infil¬ 
trated with silicon using chemical vapor deposition. The 
calculated expected widths of the bandgaps are 5 and 
7% respectively. Reflectivity was used to demonstrate 
the photonic behaviour of the crystals. Silicon inverse 
opals have been reviewed in Ref. 77|- Inverse opals of 
other high-index semiconductors such as GaAs have also 
been pursued, see Ref. [t^. 

As an alternative to opals, templates have also been 
fabricated by holographic lithography [t^ , or by interfer¬ 
ence lithography [80| where silicon inverse opals were ob¬ 
tained. Large defect channels are written in these struc¬ 
tures by an additional 2-photon polymerization process 
step. Two-photon polymerization was also used to obtain 
large waveguide-like structures in silicon inverse opals ob¬ 
tained by self-assembly [^ . 

While inverse opals are very popular on account of 
the relatively easy fabrication routes 57|, the maximum 
width of the bandgap is relatively narrow, and requires 
intricate optimization. Since the bandgap is of a higher 
order in the band structures (2nd order Bragg) [s^, the 
gap is narrower and more sensitive to unavoidable disor¬ 
der and fabrication deviations than for structures with 
a lower order band gap, which is a main disadvantage 
for inverse opals. From calculations in Ref. [s^ it is con¬ 
cluded that variations of air-sphere radii and lattice posi¬ 
tions of less than 2% of the lattice parameter are already 
sufficient to completely close the bandgap. These results 
indicate that in order to display a photonic bandgap, 
inverse opals must be made with an extremely high pre¬ 
cision, which may be beyond the present state-of-the-art 
in nanofabrication. 


Diamond-like photonic crystals 

In the early 1990s the possibility to use photonic crys¬ 
tals with diamond-like structures was described in a num¬ 
ber of papers, see Refs. [H, These photonic crystals 


are extremely interesting because they have significant 
potential for wide photonic bandgaps. Therefore, ma¬ 
terials with relatively low refractive index contrasts of 
TO > 1.9 suffice in order to obtain a photonic bandgap. 
A simultaneous advantageous feature of 3D diamond-like 
photonic crystals is that the bandgap is robust to un¬ 
avoidable fabrication deviations and random disorder, see 
Refs, [siii for calculations on this aspect. 

In one of the earliest studies calculations showed that 
air-spheres arranged in a diamond structure with a re¬ 
fractive index contrast of 3.6 would give wide bandgaps 
of up to 28% bandwidth [67|. Unfortunately such dia¬ 
mond structures remain elusive to date. In a 2004 review 
diamond-like photonic crystals and the efforts to obtain 
them have been reviewed 87[ . Here, several main results 
from that review are highlighted, as well as subsequent 
results. 



FIG. 3. Schematic illustration of a 3D woodpile photonic 
crystal that is fabricated by a layer-by-layer approach. The 
top part of the structure is lifted in order to provide a view 
on the central layer in which a point-defect with excess high- 
index material is visible. This point-defect acts as a 3D optical 
cavity, or a nanobox for light. 


Woodpiles 

3D woodpile photonic crystals were originally proposed 
by the Iowa State group as a practical way to realize 
powerful diamond-structured photonic crystal [^ . The 
expected maximum width of the bandgap was predicted 
to be a sizable 18% [1^. The crystal structure resembles 
a pile of logs of wood, hence their name. One may ar¬ 
gue that the analogy also pertains to the way they are 
often fabricated: by sequential stacking of layers of semi¬ 
conductor rods, see Fig. |3l This strategy has a distinct 
advantage: since layers are stacked in sequential fashion 
it is possible to alter the layout of individual layers. This 
freedom of design has been successfully used to incorpo- 
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rate high-quality optical cavities and waveguides in 3D 
photonic bandgap crystals, that are promising tools for 
cavity QED. In this section a few prominent examples of 
woodpiles and of woodpiles with embedded optical cavi¬ 
ties are discussed. Alternative methods to obtain wood- 
piles by direct-writing and 2-photon polymerization are 
highlighted. 

In a pioneering paper by the Sandia team, an intricate 
combination of several thin film deposition- and (silicon) 
fabrication processes was employed to obtain woodpile 
crystals [s^. Structures consisting of up to four layers 
of poly-crystalline rods were obtained. The geometri¬ 
cal properties of these woodpiles were chosen such that 
a bandgap is expected near A = 1500 nm. Measured 
spectra indicate stop gaps with a transmission reduced 
to 15% at that wavelength. While no convincing evi¬ 
dence for a bandgap was presented and the structures 
seem to have a significant misalignment between consec¬ 
utive layers, kudos are in order for this pioneering work, 
considering that a broad stopgap with a high reflectivity 
has been observed [89| . 

The Kyoto group led by Noda has been extremely suc¬ 
cessful in fabricating woodpile photonic crystals from III- 
V semiconductors and from silicon using a wafer fusion 
technique [aSIgll. In this fabrication method, photonic 
crystals are fabricated by stacking semiconductor stripes, 
see Fig. m Initially, patterns of stripes - i.e. the rods in 
the woodpile - are formed by electron-beam lithography 
and reactive ion etching in a single crystalline layer. In 
case of GaAs crystals, the layer is GaAs that sits on an 
(Al)GaAs etch-stop layer, on a GaAs substrate. A pair of 
patterned wafers is stacked and bonded in a crossed con¬ 
figuration with precise control over the alignment, and 
one of the substrates is removed. These process steps 
are repeated. At each fusing step, the number of lay¬ 
ers doubles, until the desired structure is obtained. One 
limitation of the fabrication method is that it is so ad¬ 
vanced that it has only successfully been realized by the 
Kyoto group. While initial crystals were only 4 to 8 layers 
thick, much thicker crystals were later realized. Notably, 
impressive GaAs woodpile structures have been reported 
with a thickness of 17 layers of rods [9^. To investigate 
the presence of the 3D photonic bandgap, normal inci¬ 
dence transmission and reflection spectra were collected, 
together with transmission spectra as a function of an¬ 
gle of incidence. In the latter experiments the stopband 
shifts to lower frequencies (longer wavelengths) with in¬ 
creasing incident angles, and the shift saturates at 40 to 
50°. If the range of the 3D bandgap is defined as the 
range where the attenuation exceeds 80%, it covers the 
range from 1300 to 1550 nm. 

In an elegant experiment the Kyoto group investigated 
the surface modes on a 3D bandgap crystal [o^. It was 
experimentally demonstrated that photons can be con¬ 
fined and manipulated at the surface of 3D photonic 
crystals. GaAs woodpile crystals were studied with a 
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FIG. 4. Experimental demonstration of surface states of 
a 3D photonic crystal with a 3D bandgap. (A) Experimen¬ 
tal set-up used to couple light into the surface modes. (B) 
Optical microscope image of the surface of the 3D photonic 
crystal, showing the excitation point of the incident light and 
the edge of the photonic crystal. (C-E) Experimental results 
for light propagation through the surface mode, in which the 
irradiating angle 6 was set to 45.7 deg, and the wavelength to 
Ao = 1400 nm, 1430nm, or 1460nm. (D) Light with Ao = 1430 
nm propagates throu gh t he surface modes of the 3D photonic 
bandgap. From Ref. with kind permission. 


thickness of 8 layers and a 3D bandgap at wavelengths 
between 1300 and 1600 nm. Since the surface states ap¬ 
pear in the gap at frequencies below the light line, an 
evanescent-coupling method with a prism was used to ex- 
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cite the surface states, see Fig. lUJA). Light is sent under 
total internal reflection conditions into the prism that is 
placed just above the photonic crystals. Once the evanes¬ 
cent waves excite a surface mode of the photonic crystal, 
the reflected optical power shows a distinct minimum. 
By measuring reflection spectra as a function of angle of 
incidence and tracking the reflectivity minimum versus 
wavelength, the dispersion relations of the surface states 
were obtained for various crystal directions. FigureslH^B- 
E) show images of the photonic-crystal surface observed 
from the top with a camera. Figure UfD) clearly shows 
that light with a wavelength Aq = 1430 nm propagates 
along the surface of the 3D photonic crystal, at an an¬ 
gle of incidence Q = 45.7 deg. The light propagates to 
the end of the crystal surface where it is scattered. In 
contrast, at Aq = 1400 and Aq = 1460 nm light does not 
propagate on the surface as there are no surface modes at 
these wavelengths. The wavelength of the surface modes 
could be tuned by scanning Q. This study is a creative 
demonstration of bandgap behavior, since the observa¬ 
tions reveal that the light is strictly confined to the sur¬ 
face of the crystal; this is beautiful evidence that light is 
forbidden to propagate into the bulk of the crystal. 

Ref. [ 9 ^ has also presented even more advanced ma¬ 
nipulation of ’’flat light”: by controlling the surface ter¬ 
mination, the Kyoto team was able to demonstrate gaps 
for the surface modes. Moreover, point defects were fab¬ 
ricated that behave as cavities for the surface-trapped 
light. Impressive cavity quality factors up to Q = 9000 
were reported. This work takes an important step in 
opening a new route for the manipulation of light by 
3D photonic crystals, as well as pioneering the surface 
science of 3D photonic crystals. Moreover, the results 
bear intriguing analogies to surface plasmon-polaritons 
physics [94| ; the absorption-free nature of a 3D photonic- 
crystal surface is expected to lead to new sensing appli¬ 
cations , as well as to novel light-matter interactions that 
are at the heart of cavity QED. 

In Refs. 9l|, 9^ optical cavities were realized in wood- 
pile crystals. To this end, a central layer containing point 
defects of different sizes were incorporated in the wood- 
pile stack. The intensity emitted by embedded InGaAsP 
quantum wells served to probe the cavity quality fac¬ 
tors. Several peaks were observed that could be iden¬ 
tified as cavity resonances, and a quality factor up to 
Q = 350 was measured for a thick 17-layer woodpile. 
Subsequently, GaAs and silicon woodpiles were fabri¬ 
cated with intricate waveguide structures [9i,|93. Guid¬ 
ing of light along these waveguides has been successfully 
demonstrated. These examples clearly demonstrate that 
woodpiles show great potential for control over propaga¬ 
tion and emission of light, including advanced function¬ 
ality. 

Silicon woodpiles have also been fabricated by direct 
laser writing of a template and subsequent double in¬ 
version to amorphous silicon 97|. First a woodpile was 


fabricated from photoresist, which was fully infiltrated 
with silicon dioxide through chemical vapour deposition. 
Subsequently, the resist template was removed to yield 
a silicon dioxide inverse woodpile. The filling fraction 
of the inverse woodpile can be tuned by additional de¬ 
position of silicon dioxide. In the last step, the inverse 
woodpile is infiltrated with amorphous silicon by chem¬ 
ical vapour deposition to obtain the woodpile crystal. 
Based on scanning electron micrographs, the expected 
bandgap for this structure was predicted to have a width 
of nearly 9%. Optical transmission measurements dis¬ 
played a strong stop gap in the range of the expected 
bandgap. Innovations were reported wherein the amor¬ 
phous silicon was converted to polycrystalline silicon . 
Moreover, inverse woodpiles were fabricated with embed¬ 
ded waveguides [iilgi. 

The team from Tokyo University has successfully made 
woodpiles by an impressive stacking of pre-fabricated 
layers of semiconductor rods using a micromanipulation 
technique [ 9 ^. Up to 17 layers were stacked with an ac¬ 
curacy of 50 nm. From bottom to top, the crystals consist 
of 8 GaAs layers, an active layer containing InAsSb quan¬ 
tum dots, and up to 8 GaAs layers. In some cases the rods 
in the active layer were designed such that point defects 
are incorporated, see the illustration in Figure [H For 
these cavity structures, quality factors up to Q = 2300 
were reported. It was observed that the quality factor 
increases by adding more top layers to the crystal. The 
quality factors were improved to an impressive Q = 38500 
by stacking a total o f 25 layer s and fine-tuning the size 
of the optical cavity flonllToH . 

A team from MIT has made 3D silicon photonic crys¬ 
tals that are i n essence hybrids of woodpile and air-sphere 
crystals [l02| . The structures were made with an elec¬ 
tron beam lithographic approach, where in each fabri¬ 
cated layer a hole section and rod section are vertically 
combined. E-beam lithography allows to align each sub¬ 
sequent layer to the previous one with high accuracy. 
Each layer is fabricated from deposited amorphous sil¬ 
icon and subsequent layers are deposited on top of the 
underlying ones. Since the lithographic pattern is de¬ 
fined for each layer separately, incorporation of point de¬ 
fects is possible. The calculated maximum width of the 
bandgap for these crystals is 21%. A high reflectivity up 
to 90% was observed at wavelengths around 1300 nm. 
Troughs typical of cavities resonances were observed in 
the reflectivity peaks. Unfortunately, the scanning elec¬ 
tron micrographs reveal significant fabrication-induced 
structural variations in the crystals. Nevertheless, this 
work is a beautiful example of layer-by-layer fabrication. 

It is noted that many of the woodpile photonic crys¬ 
tals discussed above appear to not be cubic and thus not 
truly diamond-like. This is caused by the fact that the 
individual stacked layers do not always possess an op¬ 
timal thickness compared to the periodicity of the rods 
in the layers. Hence the crystals have bandgap widths 
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that differ from the ones calculated for true diamond¬ 
like structures. Furthermore, layer-by-layer fabrication 
methods typically introduce relatively large alignment 
errors in the structures, which is expe cted to have an 
adverse effect on the photonic bandgap [103 1 . 


Inverse woodpiles 

Broad bandgaps with widths exceeding 25% have been 
predicted for 3D photonic crystals known as “inverse 
woodpiles” . These crystals consist of pores that run 
in two perpendicular directions in a high-index backbone. 
Thus, the structure is the inverse of the woodpile struc¬ 
ture. Compared to woodpiles, inverse woodpiles have 
several advantageous features. First, the layout and pore 
alignment is defined such that it is straightforward to 
obtain a cubic diamond-like structure with a broad gap; 
the cubic structure is not distorted by imperfect stacking. 
Secondly, the pore diameters may be varied to optimize 
the volume fraction of the high index material, which is 
an essential tuning parameter for broad bandgaps. In 
contrast, in woodpile structures an optimization of the 
volume fraction would entail a change of the nanorod 
dimensions which leads to structure distortions. 

By combining macroporous etching and focused ion 
beam milling beautiful inverse woodpiles were first fab¬ 
ricated in silicon [60 |. Impressive large structures were 
obtained although unintended misalignment reduced the 
width of the expected bandgap to 17%. These crystals 
were shown to have photonic stopbands by means of re¬ 
flectivity measurements along a high-symmetry crystal 
axis. 

Inverse woodpiles have also been fabricated by means 
of a sequence oi i) direct laser writing to obtain a poly¬ 
mer template, ii) inversion throug h dep osition of silicon, 
and in) removal of the template |104l |. The deposition 
of silicon was assisted by an intermediate silicon dioxide 
layer. The resulting structures have an expected bandgap 
width of 14% near a wavelength of 2500 nm. The occur¬ 
rence of a stop gap was confirmed by means of optical 
reflectivity and transmission measurements. 

Two-photon polymerization was us ed to obtain a tem¬ 
plate for a silicon inverse woodpile [Iq3. Inversion of 
the polymer template is achieved through a sequence of 
conformal coating of a layer of AI2O3, chemical vapor de¬ 
position of silicon, and removal of both AI 2 O 3 and tem¬ 
plate. While the demonstrated inverse woodpiles are not 
truly cubic and thus not diamond-like, a bandgap with 
a maximum width of 15% was predicted for the optimal 
structure. Measured reflectivity reveal a stop gap near 
A = 1100 nm, with maxima up to 60%. Features in the 
scanning electron micrographs and absence of interfer¬ 
ence fringes in the spectra suggest that the structures 
made by this method suffer from imperfections such as 
a nonuniform surface on large length scales and rough¬ 


ness on smaller scales, which adversely affect bandgap 
formation. 

In 2006 a two-directional etching method was intro¬ 
duced by the Kyoto grou p as a means to fabricate 3D 
inverse woodpile crystals 106|. It was shown that cubic 
diamond-like structures are obtained by reactive ion etch¬ 
ing through a suitable mask of two perpendicular sets of 
pores under angles of 45° with respect to the wafer sur¬ 
face. While the structures are thin, they show strong 
60% reflectivity peaks near 1600 nm. In 2009 a more ex¬ 
tended studied resulted in structures displaying around 
97% reflectivity [63|. In these crystals a bandgap with a 
width up to 14% is expected. In addition, clear changes 
in cw intensity were observed with embedded light emit¬ 
ters (see next section.) These results illustrate that these 
inverse woodpile are very powerful optical structures. 

A similar two-directional etching method w as u sed at 
Duke University to obtain woodpiles in GaAs 107l | . The 
maximum expected width of the bandgap of optimal 
structures is about 18%. Reflectivity experiments reveal 
a strong 92% peak near a wavelength of 1300 nm, indica¬ 
tive of a high photonic strength. An outlook was pre¬ 
sented on how to obtain microcavities at the intersection 
of unit cell modulations and line defects. 



FIG. 5. Scanning electron micrograph of a 3D inverse 
woodpile photonic bandgap crystal. Such a 3D nanostruc¬ 
ture is fabricated from monocrystalline silicon using a CMOS- 
compatible two-directional etching tech nique combined with 
advaced alingment, see Refs. [^ . Il08| . This crystal has a 
diamond-like structure. 


A new two-directional etching approach using masks 
in inclined planes was demonstrated at the University 
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of Twente to realize high quality diamond-like inverse 
woodpiles from silicon with expected bandgap widths up 
to 24% 611. The method was developed to be CMOS- 
compatible in close collaboration with high-tech indus¬ 
trial partners. In this method, high-purity single crys¬ 
talline wafers are first etched in one direction by re active 
ion etching to obtain large arrays of deep pores |109l | . Sec¬ 
ondly, the sample is cleaved, rotated by 90°, and placed 
in a dedicated holder wafer that was developed to care¬ 
fully align the samples. Thirdly, by using the holder a 
second etch mask is defined in an inclined plane to the 
first pattern with a high translational alignment accuracy 
better than 30 n m and a high rotational accuracy better 
than 0.71° |in8l |. Finally, a second set of pores is etched 
perpendicular to the first set by deep reactive ion etch¬ 
ing. The overlap region of the two perpendicular sets of 
pores form the 3D inverse woodpile crystals, see Figure[5] 
The signature of the bandgap was demonstrated by ex¬ 
tensive polarization-resolved optical reflectivity measure¬ 
ments on different crystal faces, that revealed stop bands 
that occur over large solid angles spanning I.TGtt sr [^ . 
The maximum observed reflectivity of 65% was found 
to be limited by the finite crystal thickness and by sur¬ 
face roughness. From the experimental observations, the 
bandgap was determined to have a broad relative band¬ 
width of 16%. We will see in the next section that these 
crystals are optically very powerful as they reveal promi¬ 
nent cavity QED effects on excited-state lifetimes [n3. 


other diamond-like structures 


In 1991 a type of structures with the potential for 
wide bandgaps was proposed. These structures are ob¬ 
tained by etching or milling pores of air sequentially in 
three different dir ectio ns in a backbone with a high in¬ 
dex of refraction llll . Such crystals were later fabri¬ 
cated from GaAs EH . Transmission measurements were 
used to study the optical properties of these structures 
and an attenuation up to 80% was reported. A width of 
the bandgap of around 19% was reported, inferred from 
the width of the troughs in transmittance. In addition 
this paper emphasizes the importance of avoiding making 
these structures by etching or milling of tapered pores. 
This fabrication deviation was shown to have a significant 
effect on the bandgap. 

Diamond-like silicon spiral photonic cry stals have been 
obtained by glancing angle deposition [ril| . In this 
method, silicon is grown by electron beam evaporation 
on substrates that contain tungsten seeds. These seeds 
are arranged in a suitable square lattice with a lattice 
constant near 1000 nm. Spirals were grown by carefully 
rotating the substrate during deposition. An optimal 
crystal has an expected bandgap width of nearly 15%. 
For the realized crystals a narrower bandgap width is ex¬ 
pected due to a mismatch of the obtained geometry with 


the ideal geometry. The photonic properties were ana¬ 
lyzed by optical reflectivity, where peaks up to 80% were 
observed, centered near a wavelength A = 2000 nm. 

An intriguing method to fabr icate diamond-like pho¬ 
tonic crystals is biotemplating 114| . In this method, 
the diamond-like scale of a beetle is used as a template. 
By a double, sol-gel based, inversion method, the tem¬ 
plate is replicated in titania. The periodicity of these 
diamond-like structures is in the order of the wavelength 
of visible light. The expected width of the bandgap of 
these structures is calculated to be around 5%, but the 
gap is probably reduced in width or even closed due to 
the significant structural distortions that are apparent in 
the scanning electron micrographs of the structures. 

In summary, we have seen that 3D photonic crystals 
have been fabricated by many different fabrication meth¬ 
ods. It appears that the potential and promise of wide 
bandgaps has challenged many groups all over the world 
to expand the boundaries of materials science and nan¬ 
otechnology, resulting in novel ways to sculpt and fash¬ 
ion 3D nanostructures. It is exciting to foresee how 
this know-how will be increasingly used to apply pho¬ 
tonic crystals in innovations where complete control of 
light, including single photons, is essential. For many 
of these strategies it remains a challenge to embed opti¬ 
cal cavities. In particular for strategies employing two- 
directional etching, the realization of embedded cavities 
must be demonstrated. Adapting the above - or new - 
manufacturing strategies to this end remains an inspiring 
goal for future materials science research. 



FIG. 6. Distributions of photon arrival times emitted by PbS 
quantum dots in silicon photonic crystals with a 3D photonic 
bandgap (at 0.893 eV and r = 170 nm, blue triangles), de¬ 
tuned from the bandgap (at 0.850 eV and r = 136 nm, red 
circles), and in a reference (suspension, black squares). Blue, 
black, and red curves are exponential models. The magenta 
dashed curve is calculated from the spatial distribution of the 
LDOS. The cubic crystals have lattice parameters a = 693 nm 
and c = 488 nm, wit h a = cv^, and pore radii 136 < r < 186 
nm. Data from Ref. |ll0l] . 
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CAVITY QUANTUM ELECTRODYNAMICS 

Inhibited spontaneous emission in a 3D photonic 
bandgap 

Time-resolved emission. Recently a first experiment 
was reported to study the control of the excited-stat e life - 
time of emitters in 3D photonic bandgap crystals |llO| . 
To this end, 3D inverse-woodpile crystals with a cubic 
diamond-like structure were made from silicon with fixed 
lattice parameters and a range of pore radii, so as to 
tune the bandgap relative to the spectrum of the emit¬ 
ters. As emitters PbS colloidal quantum dots were stud¬ 
ied at room temperature. The dots were immersed into 
the crystals as a dilute suspension, hence the transi¬ 
tion dipole orientations sampled all directions. The 
dots were kept in a toluene suspension to minimize non- 
radiative decay jnrad- Although toluene as a low-index 
medium reduces the refractive index contrast with silicon 
to TO = 2.3, inverse woodpile crystals have from the out¬ 
set such a broad bandgap that the gap retained a relative 
width of 5%. 

Time-correlated single photon counting was used to 
precisely measure the distribution of arrival times f{t) of 
the emitted photons [i^. Figure [S] shows time-resolved 
spontaneous emission for quantum dots in two different 
photonic crystals, compared to a reference [113. Emis¬ 
sion in a crystal outside the bandgap decays faster than 
the reference, confirming that the excited-state lifetime 
of the quantum dots is controlled by the photonic crys¬ 
tals. A limitation of this study is that the signal is not 
only emission from within the crystal, but also a back¬ 
ground of quantum dots outside the crystal. Thus at 
frequencies in the bandgap the emission is at short ar¬ 
rival times {t < 500 ns) dominated by the fast decay¬ 
ing background signal. Beyond 500 ns a slow decay is 
apparent that corresponds to a strongly inhibited emis¬ 
sion. Figure [6] reveals that in the bandgap the distribu¬ 
tion of photon arrival times decays monotonically in time 
as predicted by Ref. (s^l; no fractional decay or oscilla¬ 
tions were detected as predicted elsewhere [sl, 113, El- 
Possible reasons for this discrepancy could be that the 
predicted features are not robust to ensemble averaging 
in the experiment. Conversely, theory studies consider 
excited-state population dynamics instead of photon ar¬ 
rival times, while these two phenomena can strongly dif¬ 
fer [ 33 , or sometimes assumes an excessively large tran¬ 
sition dipole moments. 

To interpret the time-resolved spontaneous emission, 
the dynamics of the quantum dots in the crystal was 
modeled with single exponential decay, where care was 
taken to prop erly account for signal and background 
statistics jll^ - Outside the bandgap the emission rate 
was found to be up to 2-fold enhanced, in agreement 
with calculated density of states ( cf. Eq. [3|) . Within the 
photonic bandgap the emission rates are strongly inhib¬ 


ited by a factor 10 x compared to emission rates outside 
the bandgap. From additional observations it was con¬ 
cluded that the quantum dots have a very high quantum 
efficiency (> 90%), thus the low emission rates corre¬ 
spond to 10-fold enhanced excited state lifetimes, which 
is a clear step towards the stabilization of the excited 
state in a photonic bandgap. The resulting lifetime of 
Ti = 5.5^s is very long for quantum dots and good news 
for applications in quantum information processing. 

The experimental observations within the bandgap 
cannot be interpreted with current theories as these gen¬ 
erally use the plane wave expansion typical of infinite 
crystals. Therefore, a heuristic model was made for the 
spatial dependency of the LDOS in a finite photonic crys¬ 
tal, where it is postulated that the unit-cell averaged 
LDOS decreases exponentially with distance 2 : in to the 
crystal, with a new decay length called [H^ - This 

model leads to a distribution of emission rates with a min¬ 
imum rate corresponding to 160 x inhibition. The corre¬ 
sponding time-resolved emission curve is shown in Fig. [51 
where it is apparent that the curve is non-exponential due 
to the broad distribution, and the model agrees very well 
with the experiments. The estimated LDOS decay length 
^LDOS = 1.03a is 6x smaller than the Bragg length that 
can be derived from the photonic strength S (see the In¬ 
troduction chapter |ll6j| b which confirms that this is a 
new length scale typical for 3D band gaps. Finally, the 
good agreement of model and experiments suggests that 
this approach is fruitful for future ab-initio cavity QED 
theory in 3D photonic bandgaps. 

Continuous-wave intensity. In 2004 the Kyoto group 
reported the first study of emission in 3D photonic band 
gap crystals [oij. To this end, InGaAsP multiple quan¬ 
tum wells were incorporated in the central layer of 3D 
woodpile crystals with thicknesses between five and nine 
layers. Cw emission intensity spectra I{ujeg) were col¬ 
lected at room temperature as a function of the angle of 
observation. Broad stopbands were observed at constant 
angle, in agreement with broad stopgaps in the dispersion 
relation. A strongly suppressed emission was observed in 
the wavelength range between 1450 and 1600 nm, inde¬ 
pendent of the angle of observation, which agrees with the 
expected inhibition in a 3D photonic bandgap. Moreover, 
the observed range agrees well with the bandgap range 
identified by transmission measurements and by FDTD 
simulations. 

Subsequently, the Kyoto group has studied emission 
intensity in inverse woodpile crystals . A layer of In¬ 
GaAsP quantum wells was embedded in the crystal to 
serve as internal emitters. By employing advanced addi¬ 
tional wafer bonding procedures, the layer emitters was 
sandwiched between two halves of the photonic crystal, 
and carefully aligned. Emission spectra were measured 
at low temperatures (< 90 K). The cw intensity I {to eg) 
was found to reduced compared to emitters embedded in 
bulk silicon. While this observation represents a direc- 
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tional stopband, a logical next step is the demonstration 
of angle-independent emission inhibition. 

In all studies above, the reference sample was a bulk 
layer of semiconductor. Here, emitters experience a dif¬ 
ferent environment and refractive index (near 3.5) than 
in the photonic bandgap crystal (effective index prob¬ 
ably < 2). Therefore it is not possible to interpret the 
level of inhibition, since relatively large Lorentz local field 
effects might come into play Since all studies con¬ 
cern cw intensity I{uieg) whose features likely correlate 
with the LDOS, it is concluded that the emitters must 
have a relatively low quantum efficiency. This conclu¬ 
sion is borne out from the analysis of similar observations 
of inhibited cw emission in photonic crystals without a 
bandgap HP 


Intensity studies have also been performed on various 
types of emitters embedded in GaAs or Si woodpile crys¬ 
tals made by micromanipulation by the Tokyo group. Re- 
markably , no inhibited emission has been reported, see, 
e.g., Ref. [99j. It is surmised that the layer-by-layer stack¬ 
ing results in variations of the vertical layer spacings, 
which cause leaking-in of vacuum fluctuations into the 
crystal and thus adversely affect the inhibition. We con¬ 
clude that the role of fabrication imperfections on spon¬ 
taneous emission control is currently not very well known 
and therefore merits further attention. 


Emission in a nanobox for light 

An important driving force for the pursuit of 3D pho¬ 
tonic crystals is the study of emission from a cavity in 
a 3D photonic bandgap that confines light in all three 
dimensions: a nanobox for light. Pioneering step s in 

M- 


this direction were taken by the Kyoto group [91 


A central layer containing point defects of different sizes 
was incorporated in the woodpile stack. The cw inten¬ 
sity emitted by embedded InGaAsP quantum wells was 
collected to probe the optical properties of the embed¬ 
ded cavity. With increasing size of the point defect, an 
increasing number of resonances was identified, whose in¬ 
tensity has an increased baseline compared to inhibited 
emission in the pristine crystal. Apparently, increasing 
point defect size causes an increased leakage into the pho¬ 
tonic crystal and thus increasingly defies the shielding of 
the bandgap. For the smallest point defects, a single 
optical resonance was observed, whose intensity had a 
baseline equal to the signal in the bandgap range. The 
small cavities thus appear to be close realizations of a 
nanobox with a true 3D shielding by the bandgap. 

Emission in a 3D photonic bandga p cavity has also 
been studied by the Tokyo group H, llOOl . Ill7l | . In the 
first study, a layer of InAsSb quantum dots was embed¬ 
ded in the central layer where the cavity was located . 
A strong polarization-dependent variation of the cw in¬ 
tensity was reported, which could be assigned to a cavity 


resonance with a quality factor in the range of 2000. The 
main features were reproduced in a second study at low 
temperatures, w here the quality factor was boosted to 
more than 8000 |l00l |. 

In a third study, woodpile phot onic crystals with cav¬ 
ities were fabricated from silicon II7| . Ge islands were 
embedded, whose cw emission intensity was recorded at 
low temperature (25 K). Two cavity resonances with nar¬ 
row linewidths were observed. The intensity of these res¬ 
onances was reported to be 30 to 60 times greater com¬ 
pared to the reference. Since the reference consisted of 
emitters in bulk semiconductor, it is very challenging to 
extract cavity enhancement factors (Purcell factors) from 
the measurements. First, it is likely that the collection 
efficiency substantially differs between the photonic crys¬ 
tal and the bulk reference. In bulk semiconductor, the 
solid angle in which emission is collected is limited by to¬ 
tal internal reflection to a small fraction of the total solid 
angle. In contrast, light in a photonic crystal is scattered 
over all 47r solid angle, thus the collection efficiency is 
greater so that the intensity extracted from a photonic 
crystal may be overestimated. Secondly, a different ref¬ 
erence environment corresponds to a different Lorentz 
local field factor If the reference is a semiconductor, 
its refractive index is much greater, leading to a much 
higher emission rate. Hence, the intensity in the refer¬ 
ence is overestimated compared to the photonic crystal. 
Such effects require care to be properly accounted for. 

To date, studies of emission in a nanobox concern the 
cw intensity I{uieg) whose sharply peaked features likely 
correlate with the LDOS. Therefore, the emitters have a 
relatively low quan tum efficiency |37| . This conclusion is 
confirmed by Ref. [El, where it was reported that the 
Ge islands have a substantial non-radiative recombina¬ 
tion rate 'ynrad- 

We conclude that at this time, the challenge is open to 
demonstrate Purcell-enhanced emission in a 3D nanobox. 
Such a demonstration requires the use of time-resolved 
emission, and thus the embedding of quantum emitters 
with an elevated quantum efficiency into a nanobox for 
light. 


Laser action in 3D photonic crystal nanocavities 

T he p romise of miniature ’’thresholdless” laser ac¬ 
tion [20| in 3D photonic bandgap crystals has long been 
a strong motivation for the field of photonic crystals. Im¬ 
portant progress in this aspect has recently been reported 
by the Tokyo group. GaAs woodpile crystals were fab- 
ricat ed w ith a layer of point defects that act as cavi¬ 
ties 


10l| . The cavity layer was covered with an increas¬ 


ing number of crystal layers to increase the cavity con¬ 
finement to the point of achieving a quality factor up to 
38000. The cavity layer contained highly efficient InAs 
quantum dots that were excited with ns laser pulses to 
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demonstrate lasing oscillation at low temperatures. It 
was observed that the peak pump power required to reach 
the laser threshold decreased with increasing cavity qual¬ 
ity factors, as controlled by increasing the number of top 
GaAs layers. From measurements of output power ver¬ 
sus input power, the spontaneous emission coupling fac¬ 
tor was determined to be /3 = 0.54, 0.67, and 0.92 for the 
structures with 6, 8 and 12 upper layers, respectively. 
Here /3 is the fraction of spontaneous emitted light that 
contributes to lasing; it is defined as the ratio of the 
power emitted into the laser mode to the total emitted 
power over all modes. For the best confined structure, 
the value of /3 was impressively close to the theoretical 
limit of unity for a miniature thresholdless laser. In fu¬ 
ture, introducing a single quantum dot into the nanobox 
would establish an ideal solid-state system for the study 
of interactions between 3D confined photons and elec¬ 
trons enclosed in a completely controlled optical envi¬ 
ronment. 

The Tokyo group has al so f abricated woodpile pho¬ 
tonic crystals with cavities [Hi. The crystals contained 
a GaAs active layer with embedded InAs quantum dots 
that were excited with ns laser pulses to investigate las¬ 
ing at low temperatures (II K). From output power ver¬ 
sus input power measurements, the spontaneous emis¬ 
sion coupling factor was found to be as high as /3 = 0.78, 
which is probably a record for silicon microlasers. The 
coupling factor is a little lower than in the GaAs wood- 
piles mentioned above, which makes sense since the qual¬ 
ity factor is also a little lower {Q = 22000), hence the 
confinement is slightly lower. The promise for lasing in 
these photonic crystals has been clearly demonstrated in 
these impressive experiments. It is a thrilling prospect to 
see thresholdless lasing approaching room temperature, 
thereby opening prospects for the application of such in¬ 
tricate miniature lasers. 


Ultrafast all-optical switching of 3D photonic 
bandgap crystals 

It is an exciting prospect in cavity QED to rapidly 
modulate the ’’bath” that surrounds a quant um em itter, 
and thereby enter new physical regimes The 

notion to ’’switch” the density of states on ultrafast time 
scales has been first considered for 3D photonic band 
gap crystals [ 2 ^, where it was theoretically proposed to 
quickly modify the refractive index of the semiconduc¬ 
tor backbone by exciting free carriers with short laser 
pulses. As a result the 3D photonic band gap will ex¬ 
hibit a large shift in frequency and a change in width. 
At frequencies near the gap, the LDOS may be switched 
from a high value to zero, or from zero to a high value, 
or from a high to zero to a high value on 100-fs time 
scales, independent of material relaxation times. Such 
fast changes are expected to yield rich cavity QED be¬ 


havior of excited quantum emitters. To this end, rate 
equations were subsequently derived for the excited state 
population of tw o-level emitters in a time-dependent en¬ 
vironment [l2r)j | . The weak coupling approximations were 
used |12l| and the LDOS was modified on time scales 
faster than the excited-state lifetime. It was found that 
a short increase of the radiative decay rate depletes the 
excited state and drastically increases the emission inten¬ 
sity /(<) during the switch event. The time-dependent 
spontaneous emission revealed a distribution of photon 
arrival times that strongly deviated from ubiquitous ex¬ 
ponential decay: a deterministic burst of photons will 
be spontaneously emitted during the ultrashort switch 
event. 

Several experiments have been performed to study ul¬ 
trafast all-optical switching of 3D photonic bandgap crys¬ 
tals. In a pioneering study reflectivity changes we re re - 
ported on silica opaline matrices infiltrated with Si [l22l| . 
Unfortunately, this work suffered from several limita¬ 
tions: first, the low refractive index contrast was insuffi¬ 
cient for a band gap. Second, the experiments were per¬ 
formed at frequencies above the electronic band gap of 
Si, so that light is being absorbed. Finally, the maximum 
feasible refractive index change was limited by the large 
induced absorption. Soon after, the Karlsruhe gro up re - 
ported transmission changes on Si inverse opals [l23|. 
The induced absorption in their crystal was strongly re¬ 
duced by annealing the Si backbone, resulting in drasti¬ 
cally decreased Drude damping. This frequency range in 
this study was limited to the range of first order Bragg 
diffraction where a pseudogap occurs, but no pho tonic 


band gap [73|, llJ, IS^- The authors of Ref. |124| per¬ 


formed experiments on 3D Si inverse opals in the range 
of the 3D photonic ba nd g ap, and similar results were ob¬ 
tained at Minnesota 125l[ . Induced absorption was lim¬ 
ited by a judicious choice of pump conditions, allowing 
the demo nstration of a large shift of the photonic band 
gap |124l |. Fast dynamics was observed - 500 fs up and 
21 ps down - implying that switching could potentially 
be repeated at GHz rates. Ultrafast switching has also 
been performed on Si woodpile crystals that were probed 
by reflectivity over an octave in frequency including the 
telecom range [s^. Only 300 fs after the switching pulse, 
the complete band gap shifted to higher frequencies be¬ 
fore quickly relaxing within 18 ps. The switched spectra 
were successfully analyzed with a theory for finite pho¬ 
tonic crystals. 

While ultrafast switching of the optical properties of 
3D photonic crystals has been clearly demonstrated, 
the demonstration of spontaneous emission switching of 
emitters that experience a fast change of the LDOS is 
currently open. It will be an ultimate challenge to ob¬ 
serve the breaking of the weak coupling approximation 
by fast time-dependent modulation. Relevant questions 
are: what is the relevant time scale for the frequency of 
a photon in a cavity to adjust to the shifting band gap. 
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to what extent do such considerations apply to excited 
emitter states near a gap? It will be truly exciting to see 
the first light shine on this subject. 


APPLICATIONS AND PROSPECTS 

Quantum decoherence. There is currently a fast grow¬ 
ing interest in quantum information science, where the 
goal is to store, process and transmit information en¬ 
coded in inherently quantum mechanical systems 271. 
While many types of physical systems are being pursued, 
cavity QED systems involving quantum light interact¬ 
ing with quantum matter receive much attention. Solid- 
state cavity QED offers many advantageous prospects for 
qubits such as system scalability and on-chip architec¬ 
ture, miniaturization and high speed s, a nd the spatial 
localization of the quantum emitters |126l |. For the ma¬ 
nipulation of quantum states, it is paramount to pre¬ 
vent decoherence, otherwise the system will behave clas¬ 
sically 26| . It is thus desirable to increase the dephasing 
time T 2 of the system, which depends on the excited state 
lifetime Ti and the pure dephasing Tdeph within one state: 
T 2 ^ = 0. We note in passing that the in¬ 

verse of the dephasing time equals the linewidth that is 
measured in an absorption experiment. Earlier on in this 
review, we have seen that the Ti of quantum emitters is 
controlled by the LDOS, leading to already 10-fold en¬ 
hanced lifetimes observed in a 3D photonic band gap. 

Pure dephasing depends on both radiative and non- 
radiative effects. Examples of non-radiative effects are 
collisions in gas phase, vibrations in liquids, or phonons 
in solid state. These dephasing effects are typically con¬ 
trolled by cooling systems to low temperatures, or by 
mechanically decoupling the quantum system as much 
as possible from the environment by nifty spacers. The 
main source of radiative dephasing is spontaneous emis¬ 
sion from the excited state, which depends on the density 
of vacuum fluctuations, hence the LDOS. Yet again, this 
effect is amenable to control by the nanophotonic envi¬ 
ronment. In the extreme case of a photonic band gap, 
one can therefore also enhance the radiative dephasing 
time. Therefore, we conclude - in agreement with elabo- 

that a 3D pho- 


rate theoretical calculations |127l |12S 
tonic bandgap crystal offers a favorable environment to 
shield qubits operating at optical frequencies from noise 
and fluctuations. 

Resonant energy transfer. In cavity QED, it is a cen¬ 
tral question how multiple interacting emitters are con¬ 
trolled SHI- A well-known optical emitter-emitter in¬ 
teraction is resonant energy transfer between pairs of 
dipoles where an energy quantum is transferred from one 
emitter, called donor, to a second emitter, called accep¬ 
tor, see e.g. [i^. The involved dipole-dipole interactions 
are crucial to quantum information science, and (Forster) 
energy transfer plays a central role in photosynthesis, as 


well as in photovoltaics, lighting, and molecular sens¬ 
ing. Recently, the question was addressed whether energy 
transfer can be controlled by the photonic environment, 
viz. the LDOS. Dye molecules were separated by a short 
strand of double-stranded DNA, and the LDOS was con¬ 
trolled by positioning the FRET pairs near a mirror [l29l| . 
Contrary to early predictions, it was found that the en¬ 
ergy transfer efficiency does change with LDOS, whereas 
the energy transfer rate is independent of the LDOS, in 
agreement with theoretical considerations. 

It was predicted that in a 3D photonic bandgap the 
efficiency of resonance energy transfer is maximal when 
the bandgap is tuned to the donor emission frequency. 
In case of a vanishing non-radiative decay, the efficiency 
will reach a perfect 100%. The observation of changing 
energy transfer efficiency imply a change in the charac¬ 
teristic Forster distance, in contrast to common lore that 
this distance is fixed for a given pair of dipoles. In case 
of a bandgap, this interaction distance will undergo the 
greatest changes. Thus, future control of resonant en¬ 
ergy transfer in a 3D photonic bandgap promises favor¬ 
able vistas that are relevant to applications ranging from 
quantum information science to physical chemistry and 
even biophysics. 

Lighting. Photonic crystals offer the opportunity for 
spontaneous emission to be strongly controlled both in 
spatial terms (direction of emission) or in absolute terms 
(rate of emission). Therefore, it is timely to discuss 
the usage as practical light-emitting sources, including 
light emitting diodes (LED), for ev eryd ay lighting appli¬ 
cations. In a recent lucid review [l30|, it has been as¬ 
serted that photonic crystals offer favorable strategies to 
efficiently couple out the light; in comparison with usu¬ 
ally employed random surface roughening some curious 
common physical limitations were noted. The Purcell 
enhancement of the spontaneous emission rate was also 
considered as a means to enhance the internal efficiency 
of the light sources. It was found that such an approach 
is effective only for sources with specific properties, such 
as a small spatial extent and a narrow bandwidth so 
that it fits with the necessary resonant cavity. In the 
review, it is also justifiably concluded that a 3D pho¬ 
tonic bandgap is not a desired feature for lighting, as a 
bandgap results in the inevitable dec reas ing of the emis¬ 
sion quantum efficiency (see Eq. [S]) [l30|| . Interestingly, 
however, it is not widely appreciated that 3D photonic 
crystals reveal extended frequency ranges with enhanced 
density of states, which could serve to boost the emission 
rate and the (internal) quantum efficiency. Calculations 
have revealed that silicon inverse opals have 3 x enhanced 
emission rates over a broad 10% bandwidth [33|. And 
provided that one controls the orientation of the transi¬ 
tion dipole moment, a 4x enhancement is feasible over 
a huge octave-broad bandwidth (> 100% relative band¬ 
width) . Therefore, we propose that the study of 3D pho¬ 
tonic crystals at frequencies outside their bandgap may 
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yield fruitful applications. 

In the context of lighting, it is also important to con¬ 
sider the important question of how to efficiently couple 
out of a 3D photonic band gap crystal, or how to couple 
light in, e.g., to excite emitters, or to address waveg¬ 
uides and cavities. In comparison to 2D slab crystals, 
it is not very obvious how to use incoupling gratings in 
3D |l30l| . We anticipate that the powerful new method 
of optical wavefont shaping - that is, the spatial add ress- 
ing of incident waves with a spatial light modulator 131 
- will open new avenues to provide access to photonic 
crystals. In particular, based on studies where one light 
emitting particle in an optically thick scattering medium 
was selectively addressed, we are optimistic that these 
new techniques will for instance allow to effectively cou¬ 
ple light into and out of a nanobox buried deep inside a 
photonic crystal. 

Novel 3D nanofabrication strategies. In high-speed 
computing, power and especially heat dissipation issues 
are becoming increasingly important. Moreover, it is 
being realized that continued miniaturization will mee t 
boundaries set by fundamental physical limits 132l 1 133 1. 
Therefore in state-of-the-art circuitry in CMOS indus¬ 
try, an increasing number of researchers is currently con¬ 
templating to array circuits in 3D grids and networks. 
A second approach to alleviate heat loads is to perform 
data communication on even shorter length scales by op¬ 
tical instead of by electrical signals. It is our firm belief 
that 3D nanofabrication approaches that h ave been s tim- 
ulated by 3D photonic crystals fabrication will 

serve as an inspiration for novel 3D CMOS fabrication 
strategi es th at are relevant for future high-speed com¬ 
puting 132 1. 


SUMMARY 

In summary, we have seen that three-dimensional (3D) 
photonic crystals with a 3D photonic bandgap play a fun¬ 
damental role in cavity quantum electrodynamics (QED), 
especially in phenomena where the local density of opti¬ 
cal states plays a central role. One can say that pho¬ 
tonic bandgap crystals offer a knob to dial the density of 
states for broad frequency bandwidths over a wide range, 
from near zero to several times the value in vacuum. We 
have given an overview of the current status of the fab¬ 
rication of 3D photonic crystals with a bandgap at op¬ 
tical frequencies. Many different methods yield powerful 
3D photonic crystals. At this time, the widely pursued 
woodpile crystals offer the widest versatility, as embed¬ 
ded high-Q cavities and waveguides have been demon¬ 
strated. The optical experiments have been discussed 
that provide signatures of 3D bandgap behavior, such 
as broadband and wide-angle reflectivity or the obser¬ 
vation of intricate surface modes. We have discussed 
the main implications of 3D bandgaps for cavity QED, 


in particular spontaneous emission inhibition of emit¬ 
ters embedded in a 3D bandgap crystal. Here, impor¬ 
tant progress has occurred in the last decade, bringing 
inhibition from a theoretical prediction to experimental 
reality. We have discussed the progress in spontaneous 
emission and laser action of emitters placed in a photonic 
bandgap cavity. Near thresholdless laser action has been 
observed, and its realization is now approaching room 
temperature operation, thereby opening avenues for ap¬ 
plications. The steps towards the breaking of the weak- 
couping limit of cavity QED have been outlined, in par¬ 
ticular by ultrafast modulation, where experimental tools 
are steadily ripening. In the final section, we have re¬ 
viewed several exciting applications of 3D photonic band 
gaps, namely the shielding of decoherence for quantum 
optical systems, the manipulation of multiple coupled 
emitters including resonant energy transfer, lighting, and 
possible spin-off to 3D nanofabrication for future high- 
end computing. While inhibited spontaneous emission in 
a 3D photonic bandgap has tested our perseverance as 
it took twenty-five to fourty years since the predictions 
by Bykov, Yablonovitch, and John ESS, the many 
recent efforts on 3D bandgaps with favorable outcomes 
bode well for exciting contributions to nanophotonics and 
beyond. 
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TABLE II. 

Overview of various types of 3D photonic bandgap crystals that have been realized. This table provides information on the employed fabrication methods and the 
high-index backbone material of the crystals. Relevant references are provided in which detailed information on these structures and their fabrication can be found. For 
some of these references values for calculated or otherwise expected widths of the bandgaps are added between brackets. The concluding column provides additional 
remarks. The following abbreviations are used: 2PhP, two photon polymerization, CVD, chemical vapour deposition, DLW, direct-laser-writing, FIB, focused ion beam 
milling, Inv. woodpile, inverse woodpile, QD, quantum dots, RIE, reactive ion etching, SEM, scanning electron microscopy. 


Type 

Fabrication method 

Material 

References 

(indication of bandgap width) 

cQED feature 

Inverse opal 

Inversion of templates 

Si 

[7a(5%),[7^(7%),[^,[^ 

Demonstrated ultrafast 





switching 1123. 1241 

Diamond-like 





-Yablonovite 

RIE in 3 directions 

GaAs 

[1121(19%) 


-Spiral crystal 

Glancing angle deposition 

Si 

[113] (10 to 14%) 


-Biotemplated 

Double inversion of beetle scales 

Titania 

[1141(5%) 

Time-resolved emission in 





crystal with pseudogap [13^ 

-Inv. woodpile Macroporous etching -|- FIB 

Si 

[^(17%) 


-Inv. woodpile DLW -1- inversion 

Si 

ri04lfl4%l 


-Inv. woodpile 2PhP of a template -I- inversion 

Si 

[1051 f 15%) 


-Inv. woodpile 2-Directional etching 

Si 

fl06iri9%Lf63lf>14%i 

Embedded quantum well layer 

-Inv. woodpile 2-Directional etching 

GaAs 

[107] (18%) 


-Inv. woodpile 2-Directional etching with holder 

Si 

1^(24%) 

Observed inhibited and modified 





emission of QDs in 3D bandgap [11^ 

Woodpile 

Combination of many fabrication methods Si 

[^ 

Demonstrated ultrafast 





switching [8^ 

Woodpile 

Wafer fusion of layers 

III-V, Si 

[90] (16%), [9^, [92], [9^, [9^ 

Incorporated cavities and quantum wells 

Woodpile 

E-beam lithography -|- CVD 

Si 

[102] (21%) 

Incorporated point defects 

Woodpile 

DLW -I- double inversion 

Si 

[92 (8.6%), [^ (6.9%), [^ (15%), [98] 

Incorporated waveguides 

Woodpile 

Micromanipulation of layers 

GaAs, Si 

199]■ fl00hl7%L 11011.[1181.[117] f 19%) Cavity lasing observed 


to 









